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Abstract 

The Einstein static universe has played a central role in a number of emergent scenarios recently 
put forward to deal with the singular origin of the standard cosmological model. Here we study the 
existence and stability of the Einstein static solution in presence of vacuum energy corresponding 
to conformally- invariant fields. We show that the presence of vacuum energy stabilizes this solution 
by changing it to a centre equilibrium point, which is cyclically stable. This allows non-singular 
emergent cosmological models to be constructed in which initially the universe oscillates indefinitely 
about an initial Einstein static solution and is thus past eternal. 
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I. INTRODUCTION 



Recent accumulation of high resolution observations is compatible with the so called 
standard model of cosmology which has a number of intriguing features. In addition to 
the early and late accelerating phases, which are difficult to account for within the classical 
relativistic framework with non-exotic matter sources, this model also possesses an initial 
singular state at which the laws of physics break down. To deal with this latter shortcoming, 
a number of attempts have recently been made to construct models which are non-singular 



and/or past eternal 



13. 



These fall into a number of groups, including emergent scenarios 



5] (see also and cyclic/ekpyriotic models |8j]. There are also a number of 

other cosmological models which are not necessarily recurrent, but nevertheless are non- 
singular in the past. These include the pioneering non-singular model by Bojowald based on 
modifications due to Loop Quantum Gravity jj the model based on polymer matter n]] 



and others 
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The so called emergent scenarios, which are non-singular and past eternal, come in two 
varieties: those that employ classical general relativity 0, 4], and those that incorporate 
quantum effects which are expected to be present at early phases of the universe [5|. An 
important ingredient in the construction of these models is the Einstein static solution, 
which in the classical general relativistic setting has long been known to be unstable 2 . This 
instability makes the construction of emergent models within the framework of classical 
relativity difficult. Interestingly though, quantum effects which are known to be operative 
at early phases of the universe have recently been shown to be able to stabilise the Einstein 
static solution by changing it from a hyperbolic equilibrium point (which is unstable) to 



a centre equilibrium point, which is cyclically stable. This was first shown to be t 



ne case 



m presence 



also 



of quantum modifications due to Loop Quantum Gravity effects in |5| (see 



14| and other related works in this connection 15|, [l6|). In most cases, where 



the singularity is removed, the overall effect is to change the effective evolution equations 
in such a way which allow the singularity theorems to be circumvented. For example, in 



1 The history of attempts at constructing non-singular /oscillatory universes which are past eternal goes 
much further back to at least the work of Tolman [l| (see also Q and references therein for a review of 
non-singular models). 

2 This solution has, however, been shown to be stable with respect to inhomogeneous perturbations [12|). 
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the context of Loop Quantum Cosmology it has been shown (see Q| for details) that, by 

writing the effective evolution equations in terms of an effective fluid, the equation of state 

in the semi-classical regime can change, and even become phantom (with equation of state 

parameter < —1). This results in the violation of the null energy condition which removes 

a crucial barrier to singularity avoidance. 

Here we consider an ingredient not studied before, namely the contribution of vacuum 

quantum fluctuations to the energy content of classical Einstein equations, and study their 

effect on the existence and stability of the Einstein static solution. In this connection we 

□ 

recall an earlier result [17| . according to which the vacuum energy of conformal fields in a 
spatially closed static universe has a density proportional to a -4 , where a is the scale factor. 
This result has been confirmed more recently by other studies which show that the Casimir 
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We study the effects of presence of 



energy in such backgrounds is proportional to a -1 
this vacuum energy on the dynamics of the universe and show that it has the consequence 
of stabilising the Einstein static solution by changing it into a centre equilibrium point. 

The outline of the paper is as follows. In Section [III we perform a general analysis of the 
equilibrium of Einstein solution in the presence of matter and the vacuum energy. In Section 
IIHI we consider a numerical example, in which the energy content is assumed to consist of 
relativistic matter plus a vacuum term with negative pressure. In Section [IV] we study the 
existence of Einstein static solution in presence of a single scalar field, with the vacuum 
contribution represented by a self-interaction potential. We then briefly discuss possible 
alternatives for the exit from the initial Einstein phase to an inflationary one. We conclude 
with a discussion in Section [V] 

II. STUDY OF THE EINSTEIN UNIVERSE IN THE GENERAL SETTING 

In this section we study the effects of vacuum energy due to conformal fields on the 
dynamics of the universe, and in particular the way this affects the existence and the stability 
of the Einstein static solution. As was mentioned above the vacuum energy of a conformal 
scalar field in a spatially closed static universe of radius a has been shown to have a density 



given by p\ = C/a 4 , where C is a positive constant 17j |. Initially it was also proposed that 
the corresponding pressure should have the form = Pa/3, to ensure the vacuum energy- 
momentum tensor is traceless in order to respect the conformal symmetry. However, with 
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the discovery of trace anomaly such a requirement is not necessary. In the de Sitter space- 
time the equation-of-state parameter of the vacuum is uj\ = — 1, due to the symmetry of the 
background. This result, however, cannot be assumed to hold in other more general space- 
times. Here, therefore, we shall proceed by first considering the general case, p\ = ujaPa, 
where c<Ja is allowed to take arbitrary values. We shall then find the conditions for the 
stability of the Einstein static universe in terms of u\. 

Starting with a closed isotropic and homogeneous Friedmann-Lemaitre- Robertson- Walker 
model sourced by a general fluid with total density p and total pressure p, the evolution 
equations are given by 

3^ = P-|, (1) 
p + 3ff(p + p)=0, (2) 

where H is the Hubble parameter. We shall assume the fluid to consist of a combination 
of the above vacuum energy plus matter with p m = u m p m , where p m , p m and u> m are the 
corresponding density, pressure and equation-of-state parameter. Substituting these in the 
above evolution equations and letting p\ = C/a 4 , the Raychadhuri equation can be written 

as 

a = - a ^ (1 + Suj m ) + — (ui m - ui A ) ■ (3) 

The Einstein static solution is given by a = = a. To begin with we obtain the conditions 
for the existence of this solution. The scale factor in this case is given by 

2 C (U) m — Il>a) , . 

aES ~ 3^ m + l ■ (4) 

The existence condition reduces to the reality condition for a ES , which for a positive C takes 
the forms 

u m > -1/3 and uj a < u m , (5) 

or 

ijj m < —1/3 and u\ > uj m . (6) 

Therefore, in the case of ordinary matter (u m > 0) plus a positive vacuum energy with 
negative pressure (i.e. with u\ < 0), the Einstein static solution always exists 3 . 



It is interesting to note that conditions ©-([B]) exclude the case uja — 1/3, originally proposed in [17 1 
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To study the stability of this solution, it is helpful to cast the equation J3} as a 2- 
dimensional dynamical system by introducing the phase-space variables x\= a and 22 = a, 



X2 = 



Xi = x 2 , 



xl + l 



2xi 



(1 + 3u m ) + — -3 (u m - u A ) 



(7) 
(8) 



1 



In these variables the Einstein static solution corresponds to the fixed point (x\ = cles, ^2 = 
0). The stability of this equilibrium point is readily found by looking at the eigenvalues, A, 
of the Jacobian matrix = dii/dxj evaluated at this point, which are found to be 



The stability depends on the sign of A 2 . For A 2 > 0, the Einstein static solution is 
a hyperbolic fixed point and hence unstable, in the sense that trajectories starting in the 
neighbourhood of such a point exponentially diverge from it (this is the same as the classical 
relativistic case). For A 2 < 0, on the other hand, the Einstein static solution becomes a centre 
equilibrium point, which is circularly stable, in the sense that small departures from the fixed 
point results in oscillations about that point rather than exponential deviation from it. In 
this case the universe stays (oscillates) in the neighbourhood of the Einstein static solution 
indefinitely. Thus the condition for stability is given by A 2 < 0. For C > 0, this implies 
that cg>a < w m . Comparing this inequality with the conditions for existence of the Einstein 
static solution, we conclude that, given u A < u m , the Einstein universe is stable 

for u m > —1/3. In particular, it is stable in presence of ordinary matter (u m > 0) plus a 
positive vacuum energy with negative pressure 4 . 

To close this section, we note that perturbations about the fixed point imply perturbations 
in a which in turn imply perturbations in p A . In the case of vacuum energy with uo A = — 1, 
this implies a coupling between matter and vacuum. Since the conservation of total energy 
is guaranteed by (j2J), the quantum vacuum would have to exchange energy with matter. We 
shall consider this coupling in Section [IV] below, using an alternative approach. We should 
emphasise here that this coupling is essential for the above results to hold. This may explain 
why, when matter and vacuum are independently conserved and u A = —1, the stability is 

4 In the case of a negative vacuum energy, that is C < 0, the conditions for existence and stability of the 
Einstein solution are given by wa > w m > —1/3, which are satisfied by any ordinary matter, provided the 
negative vacuum term has negative pressure. 




(9) 



•ES 
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FIG. 1: The scale factor as a function of time (left) and a typical trajectory in the phase space 
(a, a) (right). 



achieved only for u m < —1/3 8|, [l5|. In the next section we shall explicitly demonstrate 



this coupling with the help of a numerical example. 



III. NUMERICAL STUDY OF THE COSMOLOGICAL DYNAMICS 

In this section we make a brief quantitative study of the effects of the vacuum energy 
on the dynamics of the universe. As an example, we consider the case where the energy 
content consists of vacuum energy, which we assume to have u>\ = — 1 and p\ = C/a 4 ", plus 
a relativistic matter with an equation-of-state parameter u rm = 1/3. Using these equation- 
of-state parameters in Eq. ([3]) we obtain 

3a 3 a + 3a 2 d 2 + 3a 2 - 2C = 0. (10) 

For the Einstein static solution the corresponding scale factor is given by a 2 = 2C/3. 
Clearly phase space trajectories starting exactly on the Einstein static fixed point would 
remain there indefinitely. On the other hand, trajectories starting in the neighbourhood 
of this point would oscillate indefinitely about this solution. In Fig. 1 we have plotted an 
example of such a universe trajectory using initial conditions given by a(0) = 1 and d(0) = 0, 
together with 2C = 3.1. 

To explicitly demonstrate the coupling between the matter and vacuum components, 
discussed above, we also plot in Fig. 2 the time dependence of the energy densities of these 
components. In the absence of coupling, we would have a constant p\ and a constant p rm a 
for relativistic matter. In our case, however, the effect of the coupling is to make these 
quantitities oscillate with time, as can be seen from Fig. 2. 
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FIG. 2: The time dependence of p rm <i 4 (left) and p\ (right). 



IV. THE SCALAR FIELD ANALOGY 

As discussed above, a perturbed Einstein universe with pa = C/a 4 and uj^ = — 1 is only- 
possible if the quantum vacuum exchanges energy with matter. In this section we consider 
an alternative procedure, originally introduced in the case of de Sitter space-time [7|], in 
which such a coupling is modelled in terms of a minimally coupled scalar field with a 
self-interaction potential V. This can be done by interpreting the usual expressions for the 
energy density and the pressure of the scalar field, 

P<t> = + V, (11) 
n = - V, (12) 

as a sum of a vacuum energy component with density V and pressure — V, plus a stiff fluid 
component with density and pressure equal to <f) 2 /2. Now, when the scale factor is perturbed, 
there is an energy exchange between the potential energy V (here representing the vacuum 
and taken to be V = C/a 4 ) and the kinetic energy 4> 2 /2. 
The Lagrangian in this case takes the usual form 

■R 1 



^-i{d4>y-v{4>) 



(13) 



For a FLRW universe with positive curvature, the evolution equations take their usual forms 

12 3 



3H 2 = V +^--^ (11) 
2 'a 2 



-H + - r (15) 



Using (ITTl) - (iT2l) . it is not difficult to cast the above system in the form of the evolution 
equations ((T])-©. Taking for the potential the expression V = C/a 4 , we can derive the 



Raychaudhuri equation in the form 



a s 'd + 2aV + 2a 2 - C = 0. (16) 

The behaviour of this system can be readily studied using the general analysis given in 
Section [Til by recalling that in this case we have a mixture of a vacuum term (with uj\ = —1) 
plus a stiff fluid (with ui m = 1). The Einstein static solution again corresponds to a = = a, 
which in this case gives the corresponding scale factor to be a 2 = C/2. A similar analysis 
to that used above shows that this fixed point is again a centre. 

It is also instructive to briefly re-visit the original classical emergent model proposed in 
j^l in terms of this alternative formulation of the dynamics in terms of a scalar field. In 
that scenario also, the initial Einstein static phase has a scalar field as its energy content. 
We can see from (|14|) - (|15p that, for the Einstein static solution, <p 2 and V are both non-zero 
constants. Therefore, while <fi is changing with time, V is not. In other words, the scalar 
field is rolling along a potential plateau. As discussed in jsj], this plateau may be considered 
as the past-asymptotic limit of a smoothly decreasing potential, which eventually leads to 
an exit from the Einstein static regime into an inflationary phase. Specific forms of such a 
potential have been considered in 

Finally, another possibility that may be considered is that of a complex scalar field. For 
example, with a harmonic field = <f)oe luJt we have cf 2 = uj 2 (pl, and equations ( |T^1) -( |T5|) are 
simultaneously satisfied, with V = w 2 ^. Therefore, V remains constant while <fi rotates in 
the complex plane. The stability of the solution indicates that this is a local minimum of 
V, and the exit to the inflationary phase may involve, for example, a tunnelling to a global 
minimum. 



V. CONCLUSION 



We have studied the existence and stability of the Einstein static universe in presence 
of vacuum energy corresponding to conformally-invariant fields. Using the result that the 
vacuum energy density in Einstein universe is proportional to the inverse fourth power of 
the scale factor, we have found the range of equation of state parameters for the vacuum 
energy such that the Einstein universe is stable, in the sense of dynamically corresponding 
to a centre equilibrium point. The importance of such a solution is due to the central role 
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it plays in the construction of non-singular emergent oscillatory models which are past 
eternal, and hence can resolve the singularity problem in the standard cosmological scenario. 

Given that the oscillatory universe discussed above is close to but not exactly an Einstein 
universe, the form of the vacuum energy density in the initial oscillatory phase of the universe 
may depart from the above inverse fourth power form. To partially answer what happens 
if vacuum density takes other forms, we considered, as a first step in this direction, a more 
general functional form of the type 

PA = C/a n . 

Proceeding in a similar manner to that used above we have been able to show that, in the 
case where the content of the universe consists of radiation (with uj m = 1/3) plus vacuum 
energy (with = —1), the Einstein static universe still exists and is stable (is a centre 
fixed point) if n > 2. 
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